Following Stellar Nucleosynthesis 








The calculation of stellar nucleosynthesis requires the simultaneous 
solution for a set of coupled differential equations, each of which has 
the form 
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To appreciate the problem, let’s examine the individual terms. 


e Term 1 represents the loss of species X in reactions with particle 
a; the value Tax (defined through (12.2.5)) is a function of temper- 
ature (through the reaction rate, ax) and particle number density 
a. The former is known through the equations of stellar structure 
(which have already been solved); an estimate of the latter is also 
known from the previous model. 


e Term 2 gives the loss of Nx due to radioactivity; specifically, 
beta-decay. The decay constant Tg is a constant of physics. 


e Term 3 behaves like term 1; it computes the creation of X due 
to nuclear reactions between b and Y. Again, the ty depends on 
temperature and the number density of b. 


e Term 4 shows the creation of X by positron decay of species Z; 
like the second term, Tp comes from atomic physics. 


Like the equations of stellar structure, the nuclear reaction equa- 
tions constitute a network of first-order, non-linear differential equa- 
tions. Unlike the stellar structure equations, however, these are not 
two-point boundary problems. Instead, they are initial value prob- 
lems, where all the abundances at t = 0 are known. Thus, the 
calculation is much more straightforward. 


Simplifications for the Solving Nuclear Reaction Network 





The entire network of nuclear reaction equations does not always 
have to be solved numerically. Many reactions quickly come to 
“equilibrium” and their rates can be expressed analytically. For 
example, consider a set of reactions where element X reacts with a 
proton to give element Y, and element Y reacts with another proton 
to yield element Z. The differential equations for these reactions 
would be 
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These could be solved simultaneously via numerical integration. 
However, suppose the reaction rate Apy is extremely high (or, equiv- 
alently, the lifetime of element Y against protons is extremely short). 
Suppose also that the initial abundance of Y was rather high. In 
this case, the derivative of (13.1.3) would be negative, and the abun- 
dance of Y would decrease rapidly. Eventually, the abundance of Y 
would fall enough so that the second term of (13.1.3) would no longer 
dominate. Equilibrium would then be reached, with dNy /dt = 0. 
When this occurs, (13.1.3) can be re-written to yield 
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Equation (13.1.3) can then be removed from the network, and Ny 
replaced by its equilibrium value in the rest of the equations. 





Nuclear Reactions in Main Sequence Stars 


The temperatures and densities in the core of a star are conducive to 
nuclear reactions. However, in main sequence stars, the reactions that 
occur do so at a rather modest rate. The reason is that the two most 
abundant species in main sequence stars (by far) are 1H and “He, yet 
collisions between these nuclei rarely result in a reaction. Specifically, 
the species 7He, °Li, and Be are all spectacularly unstable. (The 
most stable of these nuclei is Be, which has a mean lifetime of 9.7 x 
10-1” sec.) Thus 





1H + 1H —> ?He — !H+'!H 


1H + 4He —> “Li — !H + tHe 
“He + *He —> 8Be — *He + *He 


Consequently, all nucleosynthesis must occur using species which are 
relatively rare, or proceed through low-probability events. 


When one considers all the possible reactions that can occur, one 
finds that only a few are important. Most can be ignored, either 
because their cross section is very small, or because the reactants are 
extremely rare. For example 


2H + 7H —> tHe 


is a reaction that has a very high cross section, but can be ignored 
because the abundance of 7H is so low, that the probability of a single 
2H encountering another ion of 7H, rather than proton is negligible. 
If we restrict ourselves to the non-negligible reactions, then one is left 
with two main sets of reactions that change hydrogen to helium: the 
proton-proton chain, and the CNO bi- (or tri-)cycles. 


The Proton-Proton Chains 


The PP chain can be broken down into 3 subchains: 


PP-I 
1H+ 1H —> 7H+eT+y, 
7H +'H — Heg 
’He + 29He —> He + 1H+ 1H 
PP-II 
3He + “He > Be +y 
(Bese — Li+ re +y 
"Li+ H= *He + He 
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8Be — tHe + *He 
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The net result of each chain is to change two protons into two neutrons 
(and two neutrinos) and bind them up with two other protons in an 
atomic of *He. Since each hydrogen atom has an atomic mass excess 
of 7.28899 MeV, while the mass excess of 4He is 2.42475 MeV, the 
result of these chains is to liberate 26.73 MeV, or 0.7% of the rest 
mass of the proton (936.86 MeV). Most of this energy will go directly 
into the star via kinetic energy, but a few percent will be lost via the 


neutrinos. 


Of the reactions listed above, (13.2.1) is by far the slowest. It in- 
volves the weak nuclear force, and an induced positron decay of a 
proton (the rate of which must be calculation via quantum theory). 
Although it is fundamentally a different type of reaction from that 
of the strong force, its reaction rate is usually expressed in a form 
similar to the other non-resonant, strong-force nuclear reactions, t.e., 
through a value of S(E) and $2. 


Reactions involving 7H, Li, Be, and B go extremely quickly, as these 
nuclei are relatively brittle. Li, for example, can be destroyed by 
cosmic-ray spallation on the surfaces of stars, and can be burned in 
the lower envelopes of stars, far outside the core, where most nuclear 
reactions occur. Moreover, the lifetime of 8B to positron decay is 
0.78 seconds, so, for all intents and purposes, equations (13.2.8) and 
(13.2.9) can be combined into 





‘Be -+'H —+4*He+*He+e" + re +y (13.2.10) 


COMPUTING PP CHAIN NUCLEOSYNTHESIS 


The abundances of the various proton-proton chain elements are de- 


scribed by a set of differential equations. 
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where the reaction rate ri; = NiNjài j = Ni/T;(%). 


These are non-linear, coupled differential equations, which may be 
solved numerically. However, several simplifications are possible, 
which take advantage of the fact that some reactions are fast enough 
to set up equilibrium conditions. Consider the individual reaction 
rates tabulated below. 

PP Chain Reaction Rates 





Q En So da, T 
Reaction (MeV) (MeV)  KeV-barns barns (years) 
H(p, Bt v)2H 1.442 0.263 3.78x 107??? 42x107? 79x 109 
*H(p, y)?He 5.493 2.5 x 1074 7.9 x 1076 4.4 x 1078 
3He(3He, 2p)fHe 12.859 5.0 x 10° 2.4 x 10° 
3He(*He, y)/Be 1.586 4.7 x107} -28x1074 9.7x 10° 
TBe(e7, v)’ Li 0.861 0.80 3.9 x 1071 
"Li(p,4He)*He 17.347 1.2 x 10? 1.8 x 107° 
"Be(p, y)B 0.135 4.0 x 107? 6.6 x 10! 


8B(6+,v)®Be*(,4 He)*He 


18.074 7.2 3x 1078 


These numbers are from Clayton, Principles of Stellar Evolution and 
Nucleosynthesis, 1968, so are slightly out of date. The lifetimes, of 
course, are condition dependent. The values quoted here are appro- 
priate for p = 100, Tg = 15, and X = Y = 0.5. 





The first thing to notice is the extremely short lifetime of deuterium. 
Under these (solar-type) conditions, 7H will last ~ 1 hour. As a 
result, 7H quickly achieves its equilibrium abundance, 
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and N2y can be eliminated from the equations. 


Similarly “Li and “Be have typical lifetimes against hydrogen of ~ 
1 year, so they, too, will attain their equilibrium abundance. If we 
add (13.2.15) and (13.2.16), then 
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so, while 7H tracks the abundance of 1H, “Be and ‘Li will follow the 


build up of He. If we then substitute this result into (13.2.14), we are 
then left with only three coupled equations 
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Note that although these equations are more compact, they still are 
non-linear, and must be solved numerically. 


Finally, depending on the star, >He may achieve equilibrium. (The 
time it takes to do this depends sensitively on the precise reaction 
rates. At temperatures Tg < 8 K, equilibrium will take > 10° years, 
so equilibrium is usually not a good assumption in low mass stars. 
In more massive, stars, however, the time for equilibrium can be as 
short as ~ 10° years for burning temperatures of Tg ~ 20 K.) If this 
equilibrium is achieved, then setting (13.2.18) to zero results in a 
simple quadratic equation in Nsy., which has the solution 
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If you then plug this into (13.2.17) and (13.2.19), you get 








dN 

a = -cu m) Nig — 2A (3He,4He) V3 He N4He 
dN. 1 1 

~ = JACH m Nig = 5 (He, tHe) Nate sre 


or 
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dt 4 dt 


Thus, the equations are consistent. 








Energy Production in the PP Chain 


Although it is possible to simplify the PP-I, PP-II, and PP-III into 
just three (or one) equation, it is important to keep track of just 
how many reactions are occurring in each chain. Each chain helps 
turn hydrogen to helium and produces a neutrino in the process, but 
the energy carried away by the neutrino is different between the 3 
chains. The neutrino produced in the PP-I chain (13.2.1) carries 
away 0.263 MeV. Thus, if *He is synthesized entirely through PP1, 
then 2% (2 x 0.263/26.73), of the energy is lost from the system. 
Alternatively, if ?He is fused into tHe via PP-II, then a 0.80 MeV 
neutrino is produced, and 0.263 + 0.80/26.73 = > 4% of the energy 
is lost. Finally, if ?He is fused via PP-III, then a highly energetic 
7.2 MeV neutrino is created. (This is the neutrino detected in 3’Cl 
experiments.) With the PP-III chain, 27.9% of the energy is lost from 
the system. Thus, the rate of usable energy generation is 
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where F, the fraction of tHe produced via the 3 chains, is computed 
via the reaction rate ratios, t.e., 
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In equilibrium, the energy generation rate of the proton-proton chain 
is approximately 


E j 
Epp = 2:38 X 10° Y Fm Opp p X’T; 2/3 6—33.80/T5 ergs sem” 


where fpp is the electron shielding factor for hydrogen, ~ is the cor- 
rection for PPII and PPIII reactions, and 


Jpp = 1 + 0.0123 T$’? + 0.0109 T3/? + 0.0009 To 


The proton-proton chain has the least sensitivity to temperature of 
any fusion reaction. At temperatures of Tg ~ 5, v = dIn(ov)/dInT ~ 
6, while at Tẹ ~ 20, v ~ 3.5. 


Xy,:¥ = 10:90 
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Correction for neutrino energy losses due PP II and PP III chain for 
3 different compositions. 


